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ABSTRACT 

We develop a model for the regulation of galactic star formation rates Esfr in 
disk galaxies, in which ISM heating by stellar UV plays a key role. By requiring 
that thermal and (vertical) dynamical equilibrium are simultaneously satisfied 
within the diffuse gas, and that stars form at a rate proportional to the mass 
of the self-gravitating component, we obtain a prediction for Esfr as a function 
of the total gaseous surface density X and the midplane density of stars + dark 
matter p s d. The physical basis of this relationship is that the thermal pressure 
in the diffuse ISM, which is proportional to the UV heating rate and therefore to 
Esfr, must adjust until it matches the midplane pressure value set by the vertical 
gravitational field. Our model applies to regions where £ ^ 100 M pc -2 . In 
low-S SF R (outer-galaxy) regions where diffuse gas dominates, the theory predicts 
that Esfr ^ ^\^Psd- ^ ne decrease of thermal equilibrium pressure when Esfr is 
low implies, consistent with observations, that star formation can extend (with 
declining efficiency) to large radii in galaxies, rather than having a sharp cutoff 
at a fixed value of S. The main parameters entering our model are the ratio of 
thermal pressure to total pressure in the diffuse ISM, the fraction of diffuse gas 
that is in the warm phase, and the star formation timescale in self-gravitating 
clouds; all of these are (at least in principle) direct observables. At low surface 
density, our model depends on the ratio of the mean midplane FUV intensity (or 
thermal pressure in the diffuse gas) to the star formation rate, which we set based 
on Solar neighborhood values. We compare our results to recent observations, 
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showing good agreement overall for azimuthally-averaged data in a set of spiral 
galaxies. For the large flocculent spiral galaxies NGC 7331 and NGC 5055, the 
correspondence between theory and observation is remarkably close. 

Subject headings: galaxies: spiral - galaxies: ISM - ISM: kinematics and dynam- 
ics - ISM: star formation - turbulence 



Introduction 



Star formation is regulated by many physical factors, with processes from sub-pc t o 



super-kpc scales contributing to setting the overall rate (see e.g. iMcKee fc Ostrikerl 120071 ). 
One of the key factors expected to control the star formation rate is the available supply of 
gas. Over the whole range of star-forming systems, from entire spiral galaxies to circumnu- 
clear starbursts, the global average of the surface density of star formation, E SFR , is observed 
to be correlated with the global average of the neutral gas surface density E as E SFR oc E 1+p 
with 1+p 1.4 (jKennicuttll 19981 ) . Recent observations at high spatial resolution have made it 



possible to investigate local, rather than global, correlations of the star formation rate with E, 
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dent in these local studies, steeper slopes are found for the outer, atomic-dominated regions 
of spiral galaxies (as well as dwarf galaxies) compared to the inner, molecular-dominated 
regions of spirals. In addition, measured indices in the low-E, Iow-Esfr regime vary consid- 
erably from one galaxy to another. Thus, no single Schmidt law characterizes the regulation 
of star formation on local scales in the outer parts of galaxies. 

The nonlinearity of observed Schmidt laws implies that not just the quantity of gas, but 
also its physical state and the surrounding galactic environment, affect the star formation 
rate. Indices p > imply that the star formation efficiency is higher in higher-density regions, 
which are generally nearer the centers of galaxies and have shor ter dynamical times. In- 
deed , the expectation based on theory and numerical simulations (jGoldreich fc Lynden-Bell 
19651 ; iKim fc Ostrikerl l200ll ) in thin, single-phase gaseous disks is that eravi tational insta- 
bilities leading to star formation would grow only if the Toomre parameter (jToomrelll964l ) 
Q = v^k/^GH) is sufficiently small. These instabilities would develop over a timescale 
comparable to the galactic orbital time t or b = 2n/Q, which corresponds to about twice the 
two-dimensional Jeans time tj^D = Vth/{GH) when Q is near-critical. Here, Vth is the ther- 
mal speed (Vfa = -Pth/p = kT/p, for P th , p, and T the gas thermal pressure, density, and 
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temperature), and k is the epicyclic frequenc y (k 2 = R 3 dfl 2 / dR). While the implied scal- 
ing S SF r oc £f2 is roughly satisfied globally ( iKennicuttl Il998h . supporting the notion that 



Quirk! 1 19 721 ). fo r more local 



galaxies evolve towards states with Q roughly near-critical (e.g. 
observations this does not provide an accurate prediction of star formation (e.g 
20081 : 1 Wond 120091 ) . 



Leroy et al. 



In addition to galactic rotation and shear rates, an important aspect of local galactic en- 
vironment is the gravity of the stellar component. The background stellar gravity compresses 
the disk vertically (affecting the three-dimensional Jeans time tj = y/ir/ (Gp)) : and pertur- 
bations in the stellar density can act in concert with gaseous perturbations in gravitational 
instabilities (altering the effective Q). Thus, one might expect the stellar surface density 
S,,, a nd /or volume density p s to affect the star formation rate (see e.g. iKim fc Ostriker 



20071 and references therein). For example, if the stellar vertical gravity dominates that 
of the gas (see §2 for a detailed discussion of this), a scaling S SFR oc S/tj would imply 
^sfr oc S 3//2 (G 3 p s ) 1//4 /t> th for a constant-te mperature gas disk. Although star forma tion 
does appear to be correlated with E s (e.g. iRyder fc Dopital Il994j ; iHunter et al.l Il998t see 
also below), the simple scaling cx S/tj (taking into account both gaseous and stellar grav- 
ity, and assuming constant t> t h, in calcul ating tj) does not in fact provide an accurate local 
predi ction of star formation rates (e.g. 
2009h . 



Abramova k, Zasovl 120081 ; iLerov et al.l 12008c IWong 



A likely reason for the inaccuracy of the simple star formation prescriptions described 
above is that they do not account for the multiphase character of the interstellar medium 
(ISM), in which most of the volume is filled with low-density warm (or hot) gas but much 
(or even most) of the mass is found in clouds at densities two or more orders of magnitude 
greater than that of the intercloud medium. For the colder (atomic and molecular) phases, 
the turbulent velocity dispersions are much larger than v t h, so that the mean gas density p 
averaged over the disk thickness depends on the turbulent vertical velocity dispersion. Even 
when multiphase gas and turbulence (and stellar and gas gravity) are taken into account in 
simulations, however, the simple estimate S SF r oc S/tj oc T l ^/Gp (using p directly measured 
from the simulations) yields Schmidt-law indices steeper than the tr ue values measured in 
both the simulations and in real galaxies ( iKoyama fc Ostriker! l2009al ). Perhaps this should 
not be surprising, since one would expect the proportions of gas among different phases, as 
well as the overall vertical distribution, to affect the star formation rate. If, for example, 
most of the ISM's mass were in clouds of fixed internal density that formed stars at a fixed 
rate, then increasing the vertical velocity dispersion of this system of clouds would lower p 
but leave Esfr unchanged. 



The relative proportions of gas among different phases seems difficult to calculate from 
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first principles, because it depends on how self-gravitating molecular clouds form and how 
they are destroyed, both of which are very complex processes. Intriguingly, however, analysis 
of recent observations of spiral galaxies has shown that the surface density of the molecular 
component averaged over ~ kpc annuli or local patches shows a relatively simple overall 
behavior, increasing roughly linearly with the empi r ically-estimated m idplane gas pressure 
rtWong fc Blitjbooa blitz & Rosolowskylliool bood iLerov et al.lbo08h . The physical reason 
behind this empirical relation has not, however, yet been explained. 

In this paper, we use a simple physical model to analyze how the gas is partitioned into 
diffuse and self-gravitating components, based on considerations of dynamic and thermody- 
namic equilibrium. We develop the idea that the midplane pressure in the diffuse component 
must simultaneously satisfy constraints imposed by vertical force balance, and by balance 
between heating (primarily from UV) and cooling. In particular, we propose that the approx- 
im ately linear empirical relation bet ween molecular content and midplane pressure identified 
by iBlitz fc Rosolowskyi (12004 120061 ) arises because the equilibrium gas pressure is approxi- 
mately proportional to the UV heating rate; since the mean UV intensity is proportional to 
the star formation rate and the star formation rate is proportional to the molecular mass in 
normal spirals, the observed relationship naturally emerges^ We use our analysis to predict 
the dependence of the star formation rate on the local gas, stellar, and dark matter content 
of disks, and compare our predictions with observations. The analysis, including our basic 
assumptions and observational motivation for parameters that enter the theor y, is set out 
i n §2. Section 3 then compares to the observed data set previously presented in iLeroy et al. 
( 120081 ). In §4, we summarize and discuss our main results. 



2. Analysis 

2.1. Model concepts and construction 

In this section, we construct a local steady-state model for the star formation rate in the 
disk, with independent variables the total surface density of neutral gas (S), the midplane 
stellar density (p s ), and the dark matter density (pdm)- The latter two quantities enter only 
through their effect on the vertical gravitational field. To develop this model, we suppose 
that the diffuse gas filling most of the volume of the interstellar medium is in an equilibrium 
state. The equilibrium in the diffuse ISM has two aspects: force balance in the vertical 



] Dopital jl985l) previously showed that assuming the pressure is proportional to the star formation rate 
yields scaling properties similar to observed relationships. 
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direction (with a sum of pressure forces offsetting a sum of gravitational forces), and balance 
between heating and cooling (where heating is dominated by the FUV). Star- forming clouds, 
because they are self-gravitating entities at much higher pressure than their surroundings, 
are treated as separate from the space-filling diffuse ISM. The abundance of gravitationally- 
bound, star-forming clouds is nevertheless important for establishing an equilibrium state 
in the diffuse gas, because the FUV that heats the diffuse ISM originates in young OB 
associations. We assume (consistent with observations and numerical simulations) that the 
equilibrium thermal state established for the diffuse medium includes both warm and cold 
atomic gas. This hypothesis leads to a connection between the dynamical equilibrium state 
and the thermal equilibrium state: there are two separate constraints on the pressure that 
must be simultaneously satisfied. These conditions are met by an appropriate partition of 
the available neutral gas into diffuse and self-gravitating components. 

The reason for the partition between diffuse and self-gravitating gas can be understood 
by considering the physical requirements for equilibrium. The specific heating rate (r) in the 
diffuse gas is proportional to the star formation rate, which is proportional to the amount 
of gas that has settled out of the vertically-dispersed diffuse gas and collected into self- 
gravitating clouds. The specific cooling rate (nA) in the diffuse gas is proportional to the 
density and hence to the thermal pressure, which (if force balance holds) is proportional to 
the vertical gravity and to the total surface density of diffuse gas. Thus, an equilibrium state, 
in which cooling balances heating and pressure balances gravity, can be obtained by a suitable 
division of the gas mass into star-forming (gravitationally-bound) and diffuse components 
such that their ratio is proportional to the vertical gravitational field. If too large a fraction 
of the total surface density is in diffuse gas, the pressure will be too high, while the star 
formation rate will be too low. In this situation, the cooling would exceed heating, and mass 
would "drop out" of the diffuse component to produce additional star-forming gas. With 
additional star formation, the FUV intensity would raise the heating rate in the diffuse gas 
until it matches the cooling. 

In the remainder of this section, we formalize these ideas mathematically, first defining 
terms ( §2.2p . then considering the requirements of dynamical balance ( §2.3p and thermal 
balance ( §2.4p . and finally combining these to obtain an expression for the star formation 
rate when both equilibria are satisfied ( §2.5p . We then discuss, from a physical point of 
view, how the various feedback processes might act to adjust the system over time, steering 
it toward the equilibrium we have identified ( §2.6|) . While evolving to an equilibrium of this 
kind is plausible, we emphasize that this is an assumption of the present model, which must 
be tested by detailed time-dependent shnulationsg A worked example applying the model to 
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an idealized galaxy is presented in §2.71 In developing the present model, we have adopted 
a number of simplifications that a more refined treatment should address; we enumerate 
several of these issues in §2.81 



2.2. Gas components 

In this model, we divide the neutral ISM into two components. One component consists 
of the gas that is collected into gravitationally-bound clouds (GBCs) localized near the 
galactic midplane, with mean surface density (averaged over ~ kpc scales) of Egbc- The 
other component consists of gas that is diffuse (i.e. not gravitationally bound), with mean 
surface density Edifr- Here, we use the term "diffuse" in the sense of being widely dispersed or 
scattered throughout the volume; the diffuse component may include both tenuous, volume- 
filling gas and small, dense cloudlets (see below). All star formation is assumed to take place 
within the GBC component. In normal galaxies, the GBC component is identified with the 
population of giant molecular clouds (GMCs). Note that while observed GMCs in the Milky 
Way consist primarily of molecular gas, they also contain atomic gas in shielding layers. 
More generally, as we shall discuss further below, the relative proportions of molecular and 
dense atomic gas in GBCs depends on the cloud column and metallicity, and GBCs could 
even be primarily atomic if the metallicity is sufficiently low. 

The diffuse component is identified (in normal galaxies) with the atomic ISM. We treat 
the diffuse gas as a two-phase cloud-intercloud medium in thermal pressure equilibrium, with 
turbulent vertical velocity dispersion v 2 z assumed to be the same for warm and cold phases. 
Although the cold cloudlets within the diffuse component have much higher internal density 
than the warm intercloud gas, they are (by definition) each of sufficiently low mass that they 
are non-self-gravitating, such that their thermal pressure (approximately) matches that of 
their surroundings. The pressure in the i nterior of GBCs is consider ably higher than the 



pressure of the surrounding diffuse gas (cf. iKoyama fc Ostrikerll2009bf ) 



In reality, the diffuse gas would not have a single unique pressure even if the radiative 
heating rate is constant because of time-dependent dynamical effects: turbulent compres- 
sions and rarefactions heat and cool the gas, altering what would otherwise be a balance 
between radiative heating and cooling processes. Nevertheless, simulations of turbulent gas 
with atomic-ISM heating and cool ing indicate that the majority of the gas has pressure 



within ~ 50% of the mean value (IPiontek fc Ostrikerl 12005k 120071 ). although the breadth 



of the pressure peak depends on the timescale of turbulent forcing ~ -^turb/^turb compared 



Ostriker (2010, in preparation). 
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to the cooli ng time (lAudit fc Hennebellell2005[ l2010l; iHennebelle & Auditll20071 : iGazol et al. 
20051 . 120091 : Ijoung fc Mac Lowl l2006 r Goung et al.ll2009l ). Observations indicate a range of 
pressures in the cold atomic gas in the Solar neighborhood, with a small fraction of the gas 
at very high pressures, and y 50% of the gas at pressures within ~ 50% of the mean value 
rtjenkins k TripplboOll . [2007^ . 



In general, the volume-weighted mean thermal pressure at the midplane is given by 

J Pthd 3 x _ J(P th /p)pd 3 x _ J pd 3 x Jv 2 h dm 



vol 



X 



fd 3 

Po(Oi 



Id 3 



x 



J d 3 x J dm 



(1) 



where po is the volume- weighted mean midplane density of diffuse gas. The quantity (t>t h ) m ass 
is the mass-weighted mean thermal velocity dispersion; for a medium with warm and cold 
gas with respective mass fractions (in the diffuse component) f w and f c — 1 — f w and 
temperatures T w and T c , 



(0* 



fw ~\~ (1 fw) = fw 
-L in 



(2) 



Here, c w = {P w / p w ) 1 ^ 2 = (kT w / p) 1 ^ 2 is the thermal speed of warm gas. Since the ratio T w /T c 
is typically ~ 100, /„, ~ unless f w is extremely small. 

If the thermal pressure in the warm and cold diffuse gas phases are the same, (Pth)voi = 
P w = PwC w , so that from equations (|T]) and §2§, 



Pw 

Po 



(O. 



fw 



(3) 



This result still holds approximately even if the warm and cold medium pressures differ 
somewhat, since the warm gas fills most of the volume, (Pth)voi ~ Pw Note that one can 
also write p w /po = f w (Vtot/Vw) f° r Kot and V w the total and warm-medium volumes, so that 
fw ~ fw provided the warm medium fills most of the volume. If the medium is all cold 
gas, f w = T c /T w . Henceforth, we shall assume the warm and cold gas pressures are equal 
at the midplane so that (Pth)voi — > Pth', for convenience, we shall also omit the subscript on 

(' U th)mass- 



2.3. Vertical dynamical equilibrium of diffuse gas 



By averaging the momentum equation of the diffuse component horizontally and in time, 
and integrating outward from the midplane, it is straightforward to show that the difference 
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in the total vertical momentum flux across the disk thickness (i.e. between midplane and 



gdiffgnax) m ust be equal to the total weight of the d iffuse gas (e.g. iBoulares fc Coxl Il990 



Piontek fc Ostrikerl 120071 ; iKoyama &: Ostrikerl l2009bl ) . This total weight has three terms. 



The first term is the weight of the diffuse gas in its own gravitational field, 



I 



9 dz 8nG L dz ~ 2 ' 1 } 



where we have used \d^^/dz\ xemmax = 27rGS diff for a slab. The second term is the weight 
of the diffuse gas in the mean gravitational field associated with the GBCs, 



I 



Zdiff.max ^$ GBC 

P — j dz « vrGScBcSdiff, (5) 

dz 



where we have assumed that the scale height of the GBC distribution is much smaller than 
that of the diffuse gas so that \d^QBc/dz\ ~ 27tGSgbc over most of the integral. Note that 
equation (jSJ) gives an upper bound on this term in the weight, with a lower bound given by 
TrGSGBcSdiff/S, corresponding to the case in which the vertical distributions of the diffuse 
and gravitationally-bound components are the same. The third term is the weight in the 
gravitational field associated with the disk stars plus dark matter, 

2dift ' max fd$ s , d<$> dm \ , _ o . „PsdSdiff 



\ dz dz J p 

Here, p S( j = p s + Pdm is the midplane density of the stellar disk plus that of the dark 
matter halo; we have assumed a flat rotation curve V c = const, for the dark halo so that 
*. = (W/(4«G)H The stellar disk's scale height is assumed to be larger than that of 
the diffuse gas, so that g z ~ 4nGp s dZ within the diffuse gas layer. The numerical value of 
Q depends, but not sensitively, on the exact vertical distribution of the gas, which in turn 
depends on whether self- or external gravity dominates; Q ~ 0.33 within 5% for a range 
of cases between zero external gravity and zero self gravity. Allowing for a gradient in the 
vertical stellar density within the gas distribution, the stellar contribution to the weight 



3 If the vertical stellar distribution in the disk varies as p s oc sech 2 (z/H s ) with H s 



v\ s /(ttG£ s ), then the 



midplane stellar density is p s = T, S /(2H S ) — TiGTr s / (2v^ s ). Existing photometric and kinematic observations 



suggest that H s ~ const and v z 



oc 



.1/2 



(jvan der Kruit fc Searld Il982t iBottemal Il993l ) , but these are 



sensitive primarily to the central parts of the disk. Note that if the Toomre parameter for the stellar disk and 
the vertical-to-horizontal velocity dispersion ratio are both constant with radius, then p s oc p^m (V c /R) 2 . 



In this case, the ratio of the gas-to-stellar scale height is 



j/v ZlS ; since the gas can dissipate turbulence 



and cool to maintain constant while the stellar velocity dispersion secularly increases over time, the 
gas layer will tend to be thinner than the stellar layer even if both components flare in the outer parts of 
galaxies. 
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would be reduced by a factor ~ 1 — (2/3)(H g /H s ) 2 , where H g and H s are the gaseous and 
stellar scale heights. In the (unlikely) circumstance that the diffuse-gas scale height is much 
larger than that of the stars, g z ~ 2nGT, s would be substituted for the gravity of the stellar 
component, yielding a contribution analogous to that in equation (j5J) with S GB c — > 

Including both thermal and kinetic terms, and taking p — > at the top of the diffuse- 
gas layer, the difference in the gaseous vertical momentum flux between z = and Zdiff.max 
is given by P t h + Po^ 2 . The term v 2 is formally a mass- weighted quantity (analogous to 
(t>th)mass), b ut we assume a similar t urbulent velocity dispersion for the diffuse warm and cold 
atomic gas (IHeiles fc Trolandll2003l ). If the magnetic field is significa nt, a term equal to the 
difference between B 2 /(8 n) — B 2 /(4n) at z = and Zdiff,max i s added (IBoulares fc Coxl|l990 



Piontek fc Ostrikerl 20071) . Like other press ures, these magnetic terms a re volume-weighted; 
both observations ( Heiles fc Trolandl bood ) and numerical simulations f Piontek &: Ostriker 
20051 ) indicate that field strengths in the warm and cold atomic medium are similar. If the 
scale height of the magnetic field is larger than that of the diffuse gas (as some observations 
indicate; see e.g. iFerrierd 120011 ). then this term will be small, while it will provide an 
appreciable effect if B — > where p — > 0. In any case, the magnetic term in the vertical 
momentum flux may be accounted for by taking p v 2 — > p v 2 + AB 2 /2 — AB 2 = pov 2 , where 
A indicates the difference between values of the squared magnetic field at z = and Zdiff.max- 
Cosmic rays have a much larger scale height than that of the diffuse (neutral) gas, such that 
difference in the cosmic ray pressure between z = and Zdiff,max ma Y be neglected. We have 
also neglected the contribution from diffuse warm ionized gas, w hich has a low mean density 
and large scale height compared to that of the neutral gas (e.g. 



Gaensler et all 12008). 



Equating the momentum flux difference with the total weight, we have 



Pth 1 + 



c 2 f 



9 f PsdSjj 



— SdifT + ^G^GBC^diff + ^QGc^f^ 
4 -fth 



2 

diff 



(7) 



Here, we have used equations (JTJ and (j2]) to substitute fwC^/Pth for p 1 on the right-hand 
side. As noted above, the second term on the right-hand side could be reduced by up to a 
factor of two, if the scale height of the GBC distribution approaches that of the diffuse gas. 
It is convenient to define 



a = 1 



)+v 2 P th + Po v 2 + A(B 2 /2-B 2 )/(4 



7T 



fu 



Pi 



th 



which represents the midplane ratio of total effective pressure to thermal pressure. If the 
magnetic contribution is small (which would be true if AB 2 B 2 , even if magnetic and 
thermal pressures are comparable at the midplane), a is the total observed velocity dispersion 
a 2 divided by the mean thermal value. We shall treat v t , c w , and f w as parameters that 
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do not vary strongly within a galaxy or from one galaxy to another (see below), and Ediff, 
^gbc? an d Pth as (interdependent) variables. At any location in a galaxy, we shall consider 
p s d (and the total gas surface density £ = S diff + Sqbc) as a given environmental conditions. 

Equation ([7]) is a quadratic in both Sdiff and P t h- Thus, if P t h and Egbc are known, we 
may solve to obtain the surface density of diffuse gas: 

2aP th 



J diff 



GBC 



^ttGSgbc) 2 + 2wGa(P th + K d c 2 J wPsA ) 



1/2- 



(9) 



Scaling the variables to astronomical units, the result in equation 

as 



can also be expressed 



^diflf 



9.5 M pc- 2 « 



0.11 



1 M pc" 



+ 



0.011 



1 M pc~ 



+ 01 



Pth/k 
3000 Kcm- 3 



lOa/t, 



fed 



0.1 M pc" 



nl/2' 



(10) 



What are appropriate paramete r values to use? Since t hermal balance in the warm 
medium is controlled by line cooling (IWolfire et al.lll995l 120031 ). the warm medium tempera- 
ture is relatively insensitive to local conditions in a galaxy; we shall adopt c w = 8 km s _1 . Nu- 
merical simulations in multiphase gas have shown that the magnetic field is amplified by the 
magnetorotational instability to a level B 2 /(8tt) = (1 — 2)P t h, independe nt of the mass frac- 
tions of col d and warm gas and the vertical gravitational field strength (IPiontek fc Ostriker 
20051 . 120071 ) . while \B Z /B^\ <C 1. This is consistent wi th observed magnetic field strength s 
measured in the Milky Way and in external galaxies (IHeiles fc Trolandl 120051 ; iBeckl 120081 ) . 
Large-scale turbulent velocity disp ersions observed in local HI gas (both wa rm and cold) 
are ~ 7 km s _1 , comparable to c w (jHeiles &: Trolandl l2003t Mohan et al.ll2004J ). Total verti- 
cal velocity dispersions in HI gas in external galaxi es are also observed t o be in the range 
5 — 15 km s _1 , decreasing outward from the center (ITamburro et al.l 120091 ) . 



The most uncertain parameter is f w ~ f w , th e fraction of the diffuse mass in the warm 



phase. In the Solar neighborhood, this is ~ 0.6 ( IHeiles &: Trolandl 120031 ) . and in external 
galaxies the presence of both narrower and broader components of 2 1 cm emission suggests 
that both warm and cold gas are present ( Ide Blok fc Walter! 120061 ) . with some indication 
based on "universality" in li ne profile shapes that the warm-to-cold mass ratio does not 
strongly vary with position (IPetric fc Rupenl 120071 ). In the outer Milky Way, the ratio of 
HI emission to absorption appears nearly co nstant out to ~ 2 5 kpc, indicating that the 
warm-to-cold ratio does not vary significantly ( IDickey et al.ll20091 ). In dwarf gal axies as well, 
obser vations indicate that both a cold and warm HI component is present ( lYoung &: Lo 
19961 ). While uncertain, it is likely that f w ~ 0.5 — 1, at least in outer galaxies. 
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Thus, allowing for the full range of observed variation, a ~ 2 — 10, af w ~ 1 — 5, and 
a / fw ~ 2 — 20; we shall adopt a = 5 and f w = 0.5 as typical for mid-to-outer-disk conditions. 




Jenkins fc TripplhoOll and Iwolfire et allgOOgh. 



Pad. 



0.05 M fl pc- 3 ( 


Holmbere & Flvnn 


Bronfman et al.l 


1988; 


Luna et al. 



Nakanishi fc Sofud 120061 ) near the Sun, the result from equation (|T0|) is consistent 



with the obser ved total surface density estimate ~ 10 pc 2 of atomic gas in the Solar 



neighborhood ( iDickey &: Lockmanlll990l ; iKalberla fc Kerpll2009l ). 



Equation (J7J) may also be solved to obtain the thermal pressure in terms of Ediff, Sgbc, 
p s d and the diffuse-gas parameters a and f w : 



th 



4a 



Sdifi 



J GBC 
^diff 



+ 



32Crf4/ w a Psd 



ttG 



V2 
^diff 



1/2' 



'ir 



Over most of the disk in normal galaxies, the term in equation f TTTj) that is proportional to 
p s d (arising from the weight in the stellar-plus-dark matter gravitational field) dominates; 
this yields 

nCdfu 



1/2 



th 



Sditr (2Gp s d) 



1/2 



a 



(12) 



For given p s d and S, the thermal pressure therefore increases approximately proportional to 

1/2 ~ 

the fraction of gas in the diffuse phase, /die = ^dis/^- With {7rQ) ~ 1 and f w ~ f w , 
P t h ~ S«, (2Gp s d)^ 2 c w/(( v th) JrV t) 1 ^ 2 m this limit; i.e. it is the surface density of the volume- 
filling warm medium that sets the thermal pressure. Multiplying equation (fl2|) by a and 
using af w = ((v 2 h ) + v 2 )/c 2 w yields 



^difr 



(2Gp sd ) 1/2 (K 2 h )+^) 1/2 . 



(13) 



This is the same as the formula for midplane pressure adopted by iBlitz &: Rosolowskyl (120041 . 
20061 ). except that instead of Ediff their expression contains the total gas surface density S, 
instead of ((v 2 h )+v 2 Y^ 2 they use the thermal+turbulent vertical velocity dispersion (these are 
equal if vertical magnetic support is negligible), and they omit the dark matter contribution 
to p s d. Using p s = 7rGS 2 /(2f 2 8 ) and taking £ GB c Pdm — > 0, equation (ITT]) yields 



"tot ~ ~ I 1 + 



(14) 



recovering the result of lElmegreenl (119891 ) (except that he includes the total B 2 , rather than 



AB 2 , in v 2 and P t 



tot) 
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2.4. Thermal equilibrium of diffuse gas 



As expressed by equation ffTTl) . the thermal pressure in the diffuse gas must respond 
to the dynamical constraint imposed by vertical momentum conservation in the disk. In 
addition, the thermal pressure is also regulated by the microphysics of heating and cooling. 
Namely, if the atomic gas is in the two-phase regime (as is expected in a star-forming region of 
a galaxy; see §2.5p , then the thermal pressure must lie between the minimum value for which 
a cold pha se is possible, Pmin.m iri, and the maximum value for which a warm phase is possible, 



P 



max. warm ■ 



Wolfire et al.l ( 120031 ) found, based on detailed modeling of heating and cooling in 
the Solar neighborhood, that geometric mean of these two equilibrium extrema, -Ptwo-phase = 
(-Pmin,coid-Pmax,warm) 1//2 , is comparable to the local empirically-estimated thermal pressure, 
and that two phases are expected to be present in the Milky W ay out to ~ 18 kpc. Based 



on tu rbulent numerical simulations with a bistable cooling curve, iPiontek fc Ostrikerl (12005 



20071 ) found that the mean midplane pressure evolves to a value near the geometric mean 
pressure Ptwo-phase, for a wide range of vertical gravitational fields and warm-to-cold mass 
fractions. Thus, we expect the midplane thermal pressure in the diffuse gas to be comparable 
to the two-phase value defin ed by the thermal e quilibrium curve, P t h ~ Ptwo-phase- Since 
P m ax,warm/Pmin,coid ~ 2 - 5 (jWolfire et all 120031 ). even if P th = P t „o- P hase does not hold 
precisely, the midplane pressure P t h will be within a factor ~ 2 of Ptwo-phase provided the 
diffuse gas is in the two-phase regime. 



For the fiducial Solar neighborhood model of IWolfire et al.l (120031 ) . the geometric mean 
thermal pressure is Ptwo-phase A ~ Pth,oA ~ 3000 Kcm -3 . For other environments, the val- 
ues of P m in,coid and Pmax.warm depend on the heat ing of the gas: enhanced heating pushes the 
transition pressures upward (jWolfire et al.lll995l ). Because the dominant heating is provided 
by the photoelectric effect on small grains, Ptwo-phase increases approximately linearly with 
the FU V intensity. A s sumin g that Ptwo-phase scales with P m i n , C oid, we adapt the expression 



given m 



Wolfire et al.l (120031 ) and normalize Ptwo-phase using the Solar-neighborhood value: 



two— phase 
k 



12000 Kcm" 



i + 3.i(G^/a) a365 ' 



(15) 



Here, G' is the mean FUV intensity relative to the Solar neighborhood value Jfuv,o 



1 st 1 , Q is equal to the total cosmic ray/EUV/X ray ionization rate 



2.1 x 10~ 4 ergcm~ 2 s" 

relative to the value 10 -16 s -1 , and Z' d and Z' g are respectively the dust and gas abundances 
relative to Solar-neighborhood values. 

The FUV is produced by OB associations, so that the intensity will be proportional to 
the star formation rate per unit surface area. The intensity also depends on the radiative 
transfer of the UV. For example, a slab of gas with UV optical depth tj_ = kE and total 
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surface density of FUV emission Spuv has Jfuv = ^fuv(1 _ E2(t±/2))/(Attt±), where E 2 
is the second exponential integral. This correction for radiative transfer depends logarith- 
mically on 1/t\ at low optical depth (yielding a variation oc R if S is an exponential - see 
Elme green fc Parravand ( 119941 )). varying by ~ 50% for r± ~ 0.1 — 1. For simplicity, we ne- 
glect variations in Jfuv a ssociated with the radiative t ransfer here; we shall s imply assume 
Jfuv oc Spuv oc Esfr- iFuchs et al.l (120091 ) (see also iBertelli fc Nasi! 1200 ll ; IVergely et al. 
20021 ) find that the Solar-neighborhood value of the surface star formation rate Ssfr.o is 

which then yields 



2.5 x 10~ 9 M pc~ 2 yr" 1 



G' n = 



Jfuv 



J SFR 



J SFR 



" Jfuv.o ~ S SF r,o 2.5 x 10- 9 M pc" 2 yr" 1 



(16) 



Since equation (IT51) is normalized using the observed Solar-neighborhood pressure P t h,o and 
equation ( TT6l) is normalized using the (inverse of the) observed Solar-neighborhood star 
formation rate Erf r.o, our model depen ds only on the ratio of these quantities, -Pth,o/£sFR,o- 
Equivalently, since IWolfire et al.l (120031 ) predict the value of Pth,o/ Jfuv,o from theory, our 
model depends on the measured ratio of local FUV intensity to local star formation rate, 
Jfuv,o/^sfr,o- Although here we adopt an empirical value for Jfuv/^sfr based on the Solar 
neighborhood, in principle this ratio may be calculated theoretically with a detailed radiative 
transfer and population synthesis model. Simple estimates u sing standard relationships 



between the FUV emission and Esfr (e.g. ISalim et al.l 120071 ) yield values of Jfuv/^sfr 



similar to our adopted value from the Solar neighborhood. 

Assuming the high-ener gy ionization rate is proportional to the local value of Esfr and 
inversely proportional to £ (IWolfire et al.l 120031 ). 



s 



(17) 



where S is the surface density of neutral gas at the Solar circle. Strictly speaking, the above 
would only apply to cosmic rays, with £ t ' oc instead of S _1 for soft X-rays and and EUV. 
In practice, however, this does not affect the results for the star formation rate, since the 
dependence on this term is weak to begin with, and only enters the prediction for EgpR in 
outer disks (see eq. |2"2"]) . where S — > S di ff. 

Taking S ~ 10 M pc~ 2 and setting Z'JZ' = 1 yields, for P th = Ptwo-phase, 



J SFR 



6 x 10~ 10 M pc ' yr 



-2 



-1 



PthA 



3000 K cm- 



1 + 3 



10 M pc" 



o.i 



this can be combined with equation (fTTl) to yield a prediction for the star formation rate 
in terms of the gas and stellar contents of the disk. More generally, equation ( JTB"]) may be 
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inserted into equation ffT5|) . and the result substituted in equation ([9]) or f fTTj) to obtain, 
respectively, an expression for S diff in terms of S SF r, p s d, and Sgbc ; or an expression for 
E SFR in terms of S diff , p sd , and S GBC - 



2.5. The equilibrium star formation rate 



As our goal is to obtain a prediction for Esfr in terms of the total gaseous surface density 
£ = Sdiff + Sgbc and midplane stellar+dark matter density p s d, we require an additional 
relationship among the variables. Since star formation is assumed to take place only within 
GBCs, if the timescale to convert this gas to stars is ^sf,gbc; then 

^ S G bc S - S diff 

^sfr — — — • U»J 

£SF,GBC CSF.GBC 

In normal galaxies, GBCs are identified with G MCs (the outer layers of which are in fact 
atomic - see below). Recently, iBigiel et al.l (120081 ) found that there is an approximately linear 
relationship between the molecular mass measured in CO and the star formation rate, over 
the mid-disk regions in a set of spiral galaxies where X ~ 10 — 100 M pc~ 2 . The measured 
proportionality constant is £sf,co = £moi,co/£sFR ~ 2x 10 9 yr. Some recent studies targeting 
this regi me in spirals at ~ kpc scales find weak systematic variations i n tsF.Go wit h surface 



density flWong fc Blitd 12002 ; iHever et all 12004! ; iKennicutt et al.l 120071 : IVerlev et al.l |2010|). 
These are mild (no more than a factor of 2 — 3 over the range £ = 10 — 100 M pc -2 ) 
and almost all find a consistent normalization, with 2 Gyr being a typi c al timescale. CO is 
prese nt only in the portions of clouds where Ay ^ 1 (jWolfire et al.ll2010l ; iGlover fc Mac Low 
20101 ). so that it can become a poor tracer of molecular mass in low-metallicity systems. 
From observations of dust emission in the SMC, however, the star formation timescale in 
"dark" molecular gas (where hydrogen is in H2 but carbon is atomic rather than in CO) is 
found to be similar to that in CO-bright gas (A. Bolatto, personal communication). Based 
on these empirical results, £ S f,gbc = £sF,moi£ G Bc/£ m oi for ^SF.moi ~ const. 

Since GBCs contain both dense atomic shielding exteriors and dense molecular shielded 
interiors, while molecular gas may be both within GBCs and in unbound clouds, we can 
write 

SgBC -^mol,GBC + Mitom,GBC 



J mol 



Mnol,GBC + ^mol,diff 

1 + (M atom /M mol ) GB c 



(20) 



1 + S m ol,diff/S mo l,GBC 

The atomic-to-molecular ratio within externally-irradiated (spherical) clouds depends pri- 
marily on the metallicity and total cloud column density of hydrogen A^H, c ioud = M cloud / (pH^-Rcioud) 
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as ( M a tnm / M mn ]) r ] nU ri £a JZ'°- 8 (iVH, 



cloud/1-8 x 10 21 cm- 2 ) - 0.7]- 1 flKrumholz et al.l 12009a] 
McKee &: Krumholzll2010l ). Here, Nn = A^hi+2A^h2- Assuming the column densities of GBCs 
are similar to the observed values A^ H , c ioud ~ 10 22 cm -2 (i.e. S r innH ~ 100 Pc~ 2 ) typical of 



star-forming clouds in the Milky Way and Lo cal Group (e.g. ISolomon et al.lll987l ; ISheth et al. 



20081 : iBolatto et al.ll2008l : iHever et al.ll2009r ). we therefore have (M^ m /M mn ])nn c. < 1 except 



for low metallicity environments Z' < 0.2. The regions of galaxies mapped in the lBigiel et al. 
(120081 ) sample have metallicity Z' ^ 0.5, so clouds with A^ H , c ioud ~ 10 22 cm -2 would be mostly 



molecular. For clouds with this range of metallicity and column density such that Ay ^ 3, 
most of the gas that is molecular (H2 rather than H I) would also be observable in CO 
JWolfire et all I2OI0I : Idover fc Mac Low! bold ). Assuming that for moderate values of S, 
most of the molecular gas is confined within GBCs, we also have X mo i i( iifl:/£moi,GBC ^ 1 
(in galactic center regions where S > 100 M pc -2 , the diffuse molecular fraction may be 
larger). Thus, in the regime for which £sf,co is observed to be approximately constant, 
S G Bc/S mo i ~ 1 and S mol ~ X mo i,co, and we can then take £ S f,gbc ~ ^SF,moi « t SF ,co- Note 
that because CO is optically thick in clouds with A^H,cioud ~ 10 22 cm -2 and normal metal- 
licity, observed CO emission in unresolved clouds may in fact trace the atomic and "dark 
gas" portions of GBCs as well as the regions where CO is present, because these contribute 
to the gravitational potential and therefore the total CO linewidth. To the extent that this 
is true (and provided there is minimal diffuse CO-emitting gas), tsF.co would be a direct 
measurement of ^sf,gbc- 

For the rest of this paper , we shall us e tsF.n Bc —> t SF = const., and adopt the fiducial 
value of 2 x 10 9 yr based on iBigiel et al.l (120081 ). Although the value we use for ^sf,gbc 
is calibrated from observations in which star-forming clouds are primarily molecular (and 
observable in CO lines), our basic approach would remain unchanged for GBCs in different 
parameter regimes, provided that a well-defined value of ^sf,gbc is known (from either ob- 
servations with appropriate corrections for atomic and dark gas, or from theory). From a 
theoretical point of view, the internal dynamical properties of GBCs would be qualitatively 
similar whether they are mostly molecular or a mixture of cold atomic and molecular gas, but 
the value of ^sf,gbc would have to be adjusted to allow for the dependence of star formation 
efficiency on chemical content (see 



With the above assumptions, in equation §§§ we now can substitute S diff = S — Eqbc 
on the left-hand side, eliminate Pth in favor of Ssfr on the right-hand side using equation 
(Fl8|) (or eqs. [[5]-[T7]), an d set Sgbc = ^sf^sfr- The result is an implicit expression for 
SgpR in terms of S and p s d, which may be solved numerically. The Appendix provides a 
derivation of the expression for S SF r, given in equations (1A11HA12]) . Closed- form analytic 
solutions may be obtained in the limit of high and low values of the star formation rate. 
When EgpR is high, the ISM is GBC-dominated and Ssfr ~ S/isF- When Ssfr is low, the 
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ISM is diffuse-dominated. In this limit, Ediff — > X and we may drop EcBc/Sdifr in equation 
ffTTj) . resulting in 



th.low 



1700 KcnT 



10 M Q pc" 



1 + 



1 + 50f w ct' 



Psd 



0.1 M pc~ 3 



10 M pc" 2 



1/2 



(21) 



This result can then be inserted in equation ( 118]) to obtain a prediction for the dependence 
of Ssfr on p s< j and X in low-density outer-disk regions. An approximate form is 



^SFR.low 



3 x l„-° M 9 pc- yr - ( w ^) [l + 3 (^f^)" 

50/ u 



a V10 M pc 



-2 



a 



Psd 



10 M pc- 2 ^ 

l/2~ 



X 



0.1 Mq pc 



-3 



(22) 



(see also eq. IA13I in the Appendix). Note that the limit in equation ( 1221) applies in the 
Solar neighborhood, where the stellar+dark matter gravity (i.e. the term depending on p s d) 
accounts for 80% of the total contribution in second square bracket. An approximate form 
over the whole range X £ 100 M Q pc~ 2 is given by 



J SFR 



t_SF | 1 
X SsFR.low 



(23) 



(see also eq. IA14l in the Appendix). In terms of the dependence on model parameters, at high 
surface density Ssfr oc tgp, whereas at low surface density S SFR oc (/«,/«) 1//2 (EsFR,o/-Pth,o) • 



2.6. Approach to equilibrium 



The above analysis yields a relation for the star formation rate in equilibrium, but 
how would the self-consistent state of dynamical, thermal, and star formation equilibrium 
be attained in a real galaxy? First, consider the timescales involved. Vertical dynamical 
equilibrium is reached on a timescale of a few times ~ H/a z , where a z is the vertical velocity 
dispersion. With H ~ 100 pc and a z ~ 10 km s _1 , this dynamical equilibrium is typically 
reached within a few times 10 Myr. Thermal equilibrium of diffuse gas at a given density 
is accomplished by a combination of h eating and cooling; it is very fast in dense gas, and 
for low density gas requires ~ 5 Myr (IWolfire et al.ll2003l ). Separation of the diffuse gas 
into phases, which involves dynamics a nd ta kes place via thermal instability, requires ~ 
20 — 50 Myr (IPiontek fc Ostrikerl l2004j . 120051 ). The timescale for star formation to reach 
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equilibrium depends on how GBCs are formed out of the diffuse gas; recent simulations 



( IKoyama fc Ostrikerl l2009al ) suggest this may take a several tens of Myr for conditions 
similar to the Solar circle (but without spiral structure). Thus, the slowest equilibrium to be 
established is likely that of the star formation rate, which depends on the relative proportions 
of diffuse and gravitationally-bound gas. Because GBCs are continually forming and being 
destroyed within the disk (and material for a given GBC may be gathered horizontally from 
up to several H), at any time a given patch of the disk could be at a different point in this 
cycle. In an observed galaxy, measuring the equilibrium properties (for given £ and p s d) 
would require averaging over a horizontal area large enough that the different states in the 
GBC formation-destruction cycle are represented. 

As a (highly idealized) example of how the system might evolve, consider a region 
in which f^g = Sdig/S is initially high compared to the value in which all equilibrium 
conditions are satisfied. On the one hand, a higher-than-equilibrium f^is implies a lower- 
than-equilibrium proportion of gas in star-forming bound clouds, yielding a value for Esfr = 
(1 — /difi)S/tsF lower than the level when overall equilibrium (including star formation equi- 
librium) obtains. From equations (JT5l) and (|T6l) . when the radiation field is weak, the heating- 
cooling curve has a low values of -P m ax,warm and P m i njCO id (and -Ptwo-phase)- On the other hand, 
we would expect that within a few tens of Myr, thermal, dynamical, and phase equilibrium 
would be established with this (high) value of £ d iff. From equation ffTTT) . if /„, is constant, 
the midplane pressure in the diffuse gas in this situation would be higher than the value 
that would obtain when star formation equilibrium is satisfied (it can be shown algebraically 
from eq. [IT] that <9P t h/<9/difr > 0). Thus, for this situation with lower-than-equilibrium 
star- formation rate (but other processes in equilibrium), the cooling equilibrium curve sits 
at a pressure lower than that in overall equilibrium, while the midplane pressure is higher 
than that in overall equilibrium. 

Figure (PJi) shows an extreme case of higher-than-equilibrium f^m , in which the mid- 
plane pressure is higher than P m ax,warm- In this situation, a portion of the warm gas would 
condense to make cold clouds, lowering f w and reducing the midplane pressure. With a 
higher abundance of cold clouds, additional GBCs would form, lowering / diff , raising E SFR , 
and hence moving the thermal equilibrium curve toward higher pressure. As cold clouds are 
converted to GBCs, f w could return to its equilibrium value. But, with lower f diS than in the 
initial situation, the midplane pressure would be reduced. The arrows in Fig. (pji) indicate 
how the midplane pressure and thermal equilibrium curve would evolve. This process would 
continue until a self-consistent state of thermal, dynamical, and star formation equilibrium, 
with Ptwo-phase comparable to the midplane pressure P t h, is reached (see Figure [lb). (In fact, 
the midplane pressure P t h only needs to lie between Pmmxoid and P m ax,warm, which allows a 
range of self-consistent equilibria; we discuss this issue below.) 
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(a) 




P 



Fig. 1. — Schematic of approach to overall equilibrium. In (a), a higher-than-equilibrium 
value of the diffuse fraction /dis makes the midplane pressure in the diffuse gas (dashed 
line) too high, while weak heating when the bound-cloud abundance 1 — /diff is lower-than- 
equilibrium makes the thermal equilibrium curve (solid line) lie at too-low pressure. In 
( c )j /diff is lower-than-equilibrium, making the midplane pressure in the diffuse gas (dashed 
line) too low; at the same time, excess heating due to a high abundance of star-forming 
bound clouds situates the thermal equilibrium curve (solid line) too high. In (b), /diff is 
appropriate for simultaneous dynamical and thermal equilibrium with P t h at the midplane 
equal to -Ptwo-phase- Arrows in (a) and (c) indicate the direction of evolution of the midplane 
pressure level and thermal equilibrium curve for /diff to approach the equilibrium value. 
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As another example, we consider the opposite situation. We imagine a region with a 
bound cloud proportion 1 — / diff and E SFR /S above the self-consistent overall equilibrium 
values, such that heating associated with the high star formation rate makes the thermal 
equilibrium curve sit at high pressure (i.e. -Ptwo-phase is high). With low / diff , S diff = 
/diffS and the midplane diffuse-gas pressure will be low compared to the overall equilibrium 
(assuming constant f w ). Figure (JTfc) illustrates an extreme of this situation, in which the 
heating rate is so high that cold atomic clouds would secularly heat and expand to increase 
the proportion of warm gas. The energy input to bound clouds associated with the high 
star formation rate would destroy GBCs at a faster rate than they could form, increasing 
/dig and lowering Eg FR until a self-consistent equilibrium state (Figure [lb) is reached. The 
arrows indicate the evolution required move to a state in which thermal, dynamical, and star 
formation equilibrium are all satisfied. 

We note that for a given S and fdis, which fixes Sqbc and therefore Ssfr, G' and 
-Ptwo-phase for the cooling curve, the ISM conditions would still be consistent with a two-phase 
medium for all values of the midplane pressure P t h between P m i n ,coid and P max ,warm- Since 
P t h varies approximately oc (f w /a) l l 2 = f w (f w + v% /c^,)~ 1//2 (see eq. [T21) . even with f w — > 1 
the midplane pressure would increase by less than a factor of two (remaining < P max ,warm) 
compared to the solution assuming f w = 0.5. A value f w ~ 0.2 would reduce P t h by a factor 
~ 2 to approach P min , C oid- 

Although a range of f w is thermodynamically permitted, self-consistent numerical sim- 
ulations are required to ascertain how wide a range of f w can actually be realized, since 
a number of physical processes enter in setting f w . For a given /difr and £ and hence a 
given heating rate (and thus a fixed thermal equilibrium curve), solutions with a larger 
fraction of the diffuse mass in the cold phase have lower total energy; thermal instability 
and condensation of warm gas into cold clouds thus tends to drive the system toward lower 
f w . This is limited, however, by turbulent mixing and the rmal conduction, which t e nd to 



raise /,„ and also produce out-of-equilibrium gas (see e.g. iPiontek fc Ostrikerl 12005k 12007 



Audit fc Hennebelldl2005l ). Another consideration is that if the value of f w is too large (or 
too small), GBCs would not form sufficiently rapidly (or would form too rapidly) to maintain 
an equilibrium population with a total surface density S G bc — (1~ /diff)^- Yet another issue 
is that the turbulent velocity dispersion v t in the diffuse gas is likely maintained largely by 
energy inputs from star formation; with a lower star formation rate and Vt, the value of f w 
would also have to decrease in order to maintain vertical dynamical equilibrium. With self- 
consistent numerical simulations, it will be possible to assess whether f w and P t h /Ptwo-phase 
secularly depend on S and p s d- For present purposes, we proceed under the (observationally- 
motivated) assumptions that P t h w Ptwo-phase and f w ~ 0.5 — 1. The former assumption is 
also motivated by theory, as discussed in §2.41 
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2.7. Sample solution for an idealized galaxy 



As an example of the predictions from this model, we consider an idealized disk galaxy 
in which the stellar surface density obeys an exponential, T, S (R) = E s (0) exp(— R/R s ); the 
total gaseous surface density obeys a two part exponential, = S gl (0) exp(— R/R gl ) 

for the inner disk, and similarly for the outer disk with E s2 and R g 2~, and the rotation 
velocity contribution from the dark-matter halo is V(R) = V c [l — exp(—R/Rh)}. Setting 



model parameters equal to S s = 300 M pc , E s i = 150 M pc , T lg2 = 50 M Q pc , 
R s = 4 kpc, R g i = 2 kpc, R g 2 = 6 kpc, V c = 200 km s -1 , and Rh = 1 kpc, Figure (j2J) shows 
the adopted surface density profiles as well as the solution for S SFR and E diff . 



For this example, results shown in Fig. [2] (a) and 
thickness H s such that p s = S s /(0.54i?< j ) following iLeroy et al. 



c) ado p t a c onstant stellar disk 
1 The solution is 



not very sensitive to this choice, however; to show this, in Fig. [2] (b) and (d), we compare 
results for both constant-thickness stellar disk H s = const. = H g and flaring stellar disk 

_____ i i<2 

H S (R) = Hg t QE s (R s ) /T: S (R). In general, from equation (J2"2"j) . £ SFR /E oc p sd in the outer 
disk. Thus, if the stellar disk dominates, Esfr/S oc p x J 2 oc (E s /i? s ) 1//2 , which yields either 

1 /2 

Ssfr/S oc S s if H s = const., or Esfr/S oc T, s if the outer disk is flaring with v Z)S = const. 
and H s = v 2 z J (7rGS s ). If the dark matter density dominates the stellar density in the outer 
disk, then since pdm oc (V c /R) 2 , constant V c would imply E SFR /S oc R^ 1 in the outer disk. 

The example shown illust rates several chara cteris tic features that ar e in accord with the 
recent observational results of iBigiel et al.l (120081 ) and lLeroy et al.l (120081 ). First, it is evident 
that star formation does not have a sharp cutoff in the outer disk, but instead the rat e 
gradually declines with R and S s (see also the GALEX observations of lBoissier et al.l 120071 ). 
Second, there are two regimes evident for S SFR and S diff vs. S: a high-surface-density regime 
in which the gas is mostly in self-gravitating clouds and the limiting solution Ssfr = S/tsF is 
approached, and a low-surface density regime in which Ssfr has a steeper dependence. From 
equation (|22|) . the predicted limiting behavior in outer disks (or dwarfs) is S SFR oc Sp s ^ , so 
that a dropoff in the stellar and dark-matter density with R steeper than we have adopted, 
or a decline in S with R shallower than we have adopted (for a given p sd (i?)), would yield 
a steeper dependence of Ssfr on E in the outer diskj^l That is, there is no single universal 



Note that this choice of coefficient may undere stimate the midplane de nsity somewhat. For the 



local Milky Way, which has p s = 0.04 M© pc~ a (jHolmberg fc FlvnrJ l200oh and S 



42 M 



(jHolmberg fc Flvnnl 120041 ) from dynamical mass models, a scale length of R s = 2.6 kpc (jjuric et al.l 120081 ) 
would yield E s /(0.54i? s ) = 0.03 M pc" 3 . 

5 Note that in the limit of E constant, the Schmidt-law slope would be infinite, because each value of 
Psd(R) would yield a different value of Esfr- 
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Fig. 2. — Sample solution of the self-consistent star formation rate for an idealized galactic 
model: (a) shows Esfr as a function of S (solid), together with the maximum possible 
rate £/£ SF (dotted), and the outer-disk limiting solution given in equation fl22|) (dashed); 
(b) shows the ratio Ssfr/S as a function of E s , together with the maximum value 1/tsF 
(dotted); (c) shows the input stellar surface density T, s (stars) and input total gas surface 
density £ (squares) as a function of radius, together with the solution for Sdiff (triangles), 
all as a function of galactic radius R; (d) shows the solution for Esfr as a function of R. 
In (b) and (d), in addition to the case of constant stellar disk thickness (solid), results for a 
flaring disk H s oc 1/E S are shown as dashed lines. 
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slope predicted for the Schmidt law in outer disks, since p s ^ need not obey a power law 
oc TP. In individual galaxies, observations in fact show a range of behaviors for Esfr vs. E 
ranging from steep declines to slow taper at low E, with an upper envelope E SF r = E/£g F 
at high E. 

An interesting feature evident in Figure 02]) is that, even while the total gas surface 
density rises sharply in the inner disk, E di ff does not. It is straightforward to set an upper 
limit on E diff : from equation <Q, E difr < aP tiL / (-kGT, G bc)] using E G bc = E S fr£sf and 
equation ( fl8l) . 

E diff < 36 M pc- 2 Q (24) 

! + 3 (¥)' 

The physical reason for this limit is that in regions of very active star formation, the en- 
hancement in UV heating - which tends to produce diffuse atomic gas - coincides with 
stronger compression of the gas by the gravity of the massive, star-forming clouds. Under 
high-pressure conditions, the cooling rate increases; if the diffuse gas density is too high, 
cooling will exceed heating, and diffuse gas will be driven into the self-gravitating compo- 
nent. Since a ~ 2 — 10 (see eq. [H]) with a value ~ 3 — 5 most likely, the formal limit in 
equation ( |24"1) starts to become constraining only in inner galaxies where E/Eo ^> 1. As 
discussed in the Appendix (see eq. IA22I) . in fact the terms from the (variable) stellar and 
dark matter gravity in equation (J9]) reduce E diff by another factor ~ 2 relative to the upper 
limit in equation ff24|) . 

In adopting G' Q = E SFR /Eg FRi0 for equation (IT51) . we have neglected optical depth effects. 
In regions of high Edig and/or high Z' such that the optical depth is high, the mean (scaled) 
UV intensity G' within the diffuse gas is reduced (by a factor ~ [1 — E2(t±/2)]/t±; this varies 
oc 1/t_i_ for large t±), consequently lowering P th and, via equation E diff itself. Thus, the 
current simple theory overestimates E di g- in the central parts of galaxies, where tj_ becomes 
large. We further note (see §1J) that the value of E d jg is not equivalent to the surface density 
of atomic gas, since GBCs can have atomic envelopes (or even be mostly atomic at low Z) 
and since diffuse gas can be molecular (at high HZ). 



2.8. Additional considerations 

Finally, we remark on a few additional points related to assumptions behind and appli- 
cation of the theory presented above: 

1. Correction for hot gas - The formulae we have developed assume that on the scales 
over which horizontal averages are taken, the area filling factor of hot gas /hot is negligible. In 
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comparing to observations, if the fraction of the area in a beam that is primarily warm+cold 
gas is 1 — /hot, then the true surface density of warm+cold gas, for the purpose of computing 
dynamical and thermal equilibrium, is S = S obs /(l — /hot), and similarly for E diff and 
^gbc- F° r given S obs , increasing / hot would tend to increase the gaseous gravity terms in 
equation (fTTj) relative to the stellar/dark matter term, and would raise the overall value of 
the midplane pressure. 

2. The parameter f w - In the foregoing analysis, we have treated f w as an exogenous 
parameter that - based on observations - does not vary strongly within a given galaxy or 
from one galaxy to another. In reality, the relative proportions of cold and warm atomic gas 
must be self-consistently determined by considering the overall cycle of gas among phases, 
a highly complex problem. By considering a simple limiting case, however, it is possible to 
see why f w might have only moderate variation. We suppose that the primary way GBCs 
form is through self-gravitating collection of the cold diffuse atomic medium, at a rate per 
unit area £ c /t 9jC where t g>c ~ (GT Ic /(\/2ttH c ))' 1 ^ 2 is the self-gravitational timescale for the 
distribution of cold atomic clouds, which has total surface density S c « (1 — f w )Tidis and 
vertical thickness H c ~ vt/y/knGp s & (in the case where stars and dark matter dominate the 
vertical gravity) . We also suppose that GBCs have mean lifetimes £gbc (~ 20 Myr; e.g. 



Blitz et al.l 120071 ). such that the rate of destruction of the gravitationally-bound component 
per unit area is HgbcAgbc = ^sfr^sfAgbc- Equating the formation and destruction rates 
yields 

/ . v 2/3 , 2 \ 1/6 v 2/3 



Substituting in Solar-neighborhood parameters on the right-hand side, this yields f w ~ 0.6, 
in agreement with local observed estimates. Next, if we consider the diffuse-dominated 
Edig w £ limit of equation (}2"2"j) (applicability of this limit includes the Solar neighborhood) 



and substitute in Esfr oc Sy f w p s d on the right-hand side of equation ( 1251) . we see that 

(1 — f w )/ fw 3 oc (psd/S 2 ) 1 / 6 . The weak dependence on both p sd and S 2 , and the fact that 
only their ratio appears so that variations will be partially compensated, implies that f w 
would indeed be expected to vary only modestly, at least in outer disks. 

Of course, the above is an (oversimplified) description of only one of the possible ways in 
which GBCs might form. Understanding how f w depends on the fundamental environmental 
properties of a galaxy (£, p s< i, metallicity) will require numerical simulations that follow a 
wide range of processes, including realistic treatment of turbulence (which can alter f w via 
small-scale local mixing, and can also collect diffuse gas into GBCs if large-scale flows with 
long durations are present). 
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3. Galactic and metagalactic radiation - In adopting G' = £sfr/£s.f.r,o ) we have as- 
sumed that the contribution from metagalactic UV is small compared to the l ocally-generated 



i nten sity (and neglected the slowly- varying dependence on optical depth). ISternberg et al. 



(120021 ) estimate that the ratio of metagalactic to locally-generated UV in the Solar neighbor- 
hood is 0.0024. The star formation rate could therefore decline to ~ 0.002 times the local 
value, or 6xl0 -6 M kpc~ 2 yr _1 , before the metagalactic UV becomes important; this occurs 
only in the far outer regions of disks. In this regime, instead of using G' = T, S p R /T, SFR in 
equation (|T5l) . a constant value ~ 0.002 would be substituted. Beyond this point, Ptwo-phase 
would be essentially constant from equation f|T5|) . and the midplane pressure (given by eq. [2~T1) 
would fall below P m i n , C oid (as in Fig. [lb). Beyond this point, the diffuse gas would be essen- 
tially all warm, and star formation could only occur to the extent that gas is externally com- 
pressed, e.g. by spiral density waves. In practice, since in our theory S SFR oc E^/psd cx HV C /R 
in outer disks if dark matter domi nates the gravity, the gas may reach sufficiently low sur- 



face density (~ 0.1 M pc 2 ; see ISternberg et al.l 120021 ) that it could be ionized by meta 



galactic X-rays before the point where £ ;S 0.002 M Q pc~ 2 R/ kpc is reached. In addition 
to true metagalactic radiation, the contribution to the UV intensity originating nonlocally 
within a galaxy climbs as the local optical depth drops. At very low optical depths where 
Jfuv oc SsFRln(l/rj_), the star formation rate would be reduced below that in equation (1221) 
by a factor ~ 1/ ln(l/r^); this varies oc 1/R if the gas surface density obeys an exponential. 
In far outer galaxy regions, the radiation from the inner galaxy may exceed that produced 
by local star formation, expecially if the disk is strongly flaring. 

4. Properties of GBCs - In this work, we have not focused on the details of individual 
GBCs, beyond making the simplifying assumption that the gas in these structures forms 
stars on a timescale tsF that does not sign i ficant ly vary from one cloud to another. At 



a fundamental level (see iMcKee fc Ostrikerl 120071 ). the star formation timescale within a 



GBC is expected to depend on the mean density (which sets the mean gravitational free-fall 
time tfj) and on the amplitude of turbulence and strength of the magnetic field (since these 
properties determine how gas is further compressed and rarefied). 

If the star formation rate in a GBC is defined to be £ffM cloud /t ff , then £sf,gbc —ts/ £$, 
wher e t ff = (ttMh^h/SL. . J 1/4 ( 8G)~ 1/2 in terms of the cloud's mass and mean surface den- 



sity. iKrumholz fc McKed (120051 ) have argued that, due to the lognormal form of the density 
PDF in turbulent clouds, Eg will depend only weakly on a cloud's internal turbulent Mach 
number, which itself varies as fturb/^th oc (Mdo^Tido^y^T^^. As a consquence, t^^nc, 



i s not expected to vary very strongly with a cloud's properties; e.g. IKrumholz fc McKee 
( 2005 ) propose a scaling which yields tsF.GBC oc M cl ^ d £ cl 2 ^ T~jj^. As noted above, GBCs 



are composed of both molecular (shielded) and atomic (shielding) gas. Because the atomic 
gas has temperature somewhat higher ( 100 K) than that of the molecular gas ( ;$ 10 K), 
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the star formation efficiency may vary within a given GBC, as well as varying from one 
cloud to another. Since the temperature is determined by cooling, it depends on whether 
carbon is mostly atomic o r in CO, with the latter holding in the more-shielded parts of 
clouds ( IWolfire et al.ll2010l ). In a more refined theory, these intra- and inter-cloud variations 
could be taken into account in determining a mean value of £sf,gbc (for an assumed cloud 
mass function); here, we have simply adopted a single constant value, tsF- 

It is worth emphasizing again that the GBC component in our model is not equivalent 
to the molecular component observed in galaxies. An individual GBC is composed of a 
mixture of molecular gas and cold atomic gas that depends on shielding, and could be 
primarily atomic at sufficiently low met allicity. For spherical c louds, {M^m I M mn \) r ] m( \ & 
[Z'°- 8 (N H<cloud /1.8 x 10 21 cm" 2 ) -0.7]- 1 fcrumholz et alil2009al : TMcKee fc Krumholzlboioh . 



Although the chemistry of GBCs depends strongly on Z, the temperature o f the cold gas 



(e.g. 



Wolfire et al. 



is rel atively insensitive to metallicity at high density Uh ~ 100 cm 
19951 ). Thus, the internal dynamics of primarily-atomic GBCs - including the processes 
that determine the internal star formation efficiency - are expected to be similar to those in 
primarily-molecular GBCs, provided that their gravitational potentials and internal velocity 
dispersions are similar so that t> tur b fth,coid- 



We note that iKrumholz et al.l (l2009bf ) have developed a model for galactic star for- 
mation rates under an alternative set of assumptions. In their model, interstellar gas is 
assumed to be gathered into complexes with mean surface densities S comp i ex ~ 5S, where 
£ is averaged over ~ kpc scales. The fraction of mass within complexes that participates 
in star formation is determined by shielding. The molecular gas is assumed to be in GMCs 
with a surface density equal to the value observed in local galaxies, S c ioud = 85 M© pc~ 2 , 
provided this excee ds the mean gas sur f ace de nsity. Stars form in the GMCs at a rate 
determined by the IKrumholz fc McKed (120051 ) theory. In this model, Ssfr becomes a 
steep function of £ when complexes be come primarily a tomic, for ISM surface density 
£ ^ (20/Z')(£/£ comp i cx ) M pc~ 2 . The Krumholz et al. (j2009b ) model makes fewer as- 
sumptions and extends to higher-£ conditions than the model discussed here. However, 
while it is successful in describing the average properties of star-forming galaxies, it is sub- 
stantially less accurate than the present model in describing the star formation in individual 
galaxies, which is discussed below. 



3. Comparison to Observations 



The formulae derived above yield predictions for Esfr as a function of galactic gas and 
stellar properties, and can be compared with observations. Here, we compare with a recent 
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survey of s piral galaxies , for w hich the gas, stellar, and star formation content is described 



in detail in iLeroy et al.l (120081 ) . The observed measurements include molecular and atomic 
gaseous surface densities, E mo i and S atom (based on CO J = 2 — 1 and 1 — maps and 
on HI 21 cm maps, respectively, and corrected for Helium), rotation curves V C (R), stellar 
surface densities E s (based on 3.6 pm Spitzer maps), and star formation surface densities 
Esfr (based on FUV maps from GALEX and 24 /im Spitzer maps). The surveyed regions 
include both molecule-dominated and atomic-dominated areas, extending to ~ 1.2r25- 

Here, we estimate stellar densities in t wo ways, taking th e disk scale height H s = const. 



so that p s (R) = £ s (-R)/(0.54-R s ) following iLerov et al.l (120081 ); and taking H s oc 1/E S (i.e. a 



flared disk) so that p s (R) = T%(R)/[0.54:R s T ls (R s )], where R s is the fitted exponential scale 
length of the stellar disk. We estimate dark matter densities using observed rotation curves, 
as Pdm(-R) = (Yc ~ KfsV^^G-R 2 ), where V^ s /R is the correction for the contribution to the 
radial acceleration from the stellar disk (this correction is ^50% in the outer disk, where 
the contribution to vertical gravity from dark matter becomes significant). Given p s , p^m, 
and total gas surface density E = E atom + E mo i at each radius R, we numerically solve the 
equations developed in §2 to obtain predictions for Ediff and Esfr (see the Appendix for 
details). For the initial comparisons presented here, we use annular averages of the data sets 
in each galaxy; "pointwise" comparisons we have made at map resolutions of 800 pc show 
similar results, in terms of the mean values and the scatter in the observations and model 
predictions. 

For the galaxy NGC 7331, Figures [3] and [4] present results for model predictions com- 
pared to the observations, for flat and flared stellar disk cases H s = const, an d H x oc !/£,«,, 



respec tively. For this galaxy, we have used the fitted metallicity profile of iDutil fc Roy 
(119991 ) and assumed Z' d = Z' , for equation (fTBj) . We have adopted the fiducial values a = 5 



and f w = 0.5 for the diffuse-ISM parameters, as discussed in §2. The values of Esfr and 
Esfr/S are shown as functions of R, E, and E s . Also shown are the input profiles of E s , 
E = E atom + E mo i, and E atom , together with the predicted E difI . Evidently, the model pro- 
vides a remarkably good prediction for Esfr, with a slightly better match for the flared-disk 
case ( £ 20% differences) than for the flat-disk case ( £ 50% differences). In particular, the 
prediction follows the observation quite well in the atomic-dominated regime (see Fig. Hj), 
outside of ~ 7 kpc. We note that with slight adjustments of f w /a, the flaring of the stellar 
disk, or the dark matter density compared to the standard parameters and prescriptions, 
even closer agreement between predicted and observed Esfr can be obtained. Figures [3] 
and H] show that the model result for E^g exceeds Ehi in the inner region of the galaxy. 
While some of the diffuse gas in galactic center regions could in fact be molecular, we note 
(see §2.7p that our neglect of radiative transfer in estimating Jfuv makes Edig increasingly 
inaccurate in regions of high E. This does not affect the predicted value of Esfr, however, 
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Fig. 3. — Comparison between observed annular averages and model prediction for the galaxy 
NGC 7331. In (a), (b), and (d), filled squares with error bars show the observed values of 
SgFR or Eg FR /E, where E is the total gas surface density; open triangles (red) show the 
predictions from the theory of §2. In (c), stars show the observed stellar surface density E s , 
squares show total observed E, filled triangles show the observed atomic gas surface density, 
and open triangles (red) show the model prediction for the diffuse-gas surface density E diff . 
For the observations, error bars indicate scatter in the values within each azimuthal ring 
(systematic errors are larger). For the model prediction, H s = const, is adopted for the 
stellar scale height. 
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Fig. 4. — Same as in Fig. EJ except a flaring stellar disk with H s oc 1/E 8 is adopted. 
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because Edig <C S. 

For the remaining set of ten spiral galaxies in the sample, Figures and [6] show the 
comparison between observed and predicted values of S SFR , for H s = const, and H s oc 1/E S . 
As for NGC 7331, we adopt the fiducial parameter values f w = 0.5 and a = 5; no corrections 
for departures from Solar gas-to-dust ratio have been made, however. Overall, the predictions 
follow the observed profiles fairly well. For NGC 5055, which, like NGC 7331, is a flocculent 
galaxy, quite good agreement is evident ( £ 30% differences for H s = const, or £ 50% 
differences for H s oc 1/S S ). Again, this is true both for the inner molecule-dominated 
region, and for the atomic-dominated region which lies outside of ~ 8 kpc. When considering 
azimuthally-averaged data, it is not surprising that the prediction should be most accurate for 
flocculent galaxies: since the dependence of Ssfr on £ is nonlinear at high E, the prediction 
for the average Xsfr using the azimuthally-averaged S can significantly underestimate the 
average using local values of £ if the gas and star formation are both strongly concentrated 
in spiral arms (see eq. IA20l in the Appendix). Note that both NGC 7331 and NGC 5055 are 
also observed to have particularly clean linear relationships between Esfr and £ mo i. 

For a few cases (e.g. NGC 3198, NGC 4736, NGC 5194), the shape of the predicted S SFR 
follows that of the observations in the atomic-dominated part, but with an offset downward 
in the overall magnitude. From equation f l22|) . Ssfr oc (fw/ct) 1 / 2 in the diffuse-gas dominated 
regime, so that adopting a number for f w /a larger than the fiducial value we have assumed 
would shift the predicted S SF r upward. Of course, it is also possible that observational 
systematics contribute to this offset. For NGC 4736, the inner molecular-dominated regime 
has an offset between predicted and observed log(EsFR) that is similar to the outer-disk 
offset. Since the predicted Esfr in inner disks is oc t§p but independent of f w and a, the 
similarity of the inner- and outer-disk offsets could either mean that tgp and (f w / a) 1//2 happen 
to vary together, or that the observed star formation rate is systematically overestimated 
everywhere (n ote that this galaxy has other peculiarities in both its star formation and gas 



dynamics; see lWong fc Blitzll2000l ). The galaxies NGC 3198 and 3351 have observed values 
of S/Ssfr in their inner regions less than the fiducial value tsF = 2 x 10 9 yr; potentially, the 
gas surface density may be underestimated if the CO-to-H 2 conversion factor is too small, 
or star formation rates may be overestimated. 

For some other galaxies (e.g. NGC 0628, NGC 6946, and to a lesser extent, NGC 3184) 
the shape of the predicted and observed Ssfr differ somewhat. The sense of the discrepancy 
is that the observed Esfr is higher than the predicted value at intermediate radii. If some 
star- forming gas is present that is not observable either in 21 cm or CO lines, this could in 
part account for the discrepancy. These galaxies also have an irregular - and sublinear on 
average - relationship between Ssfr and £ mo i as inferred from CO. This might indicate that 
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Fig. 5. — Comparison between observed annular averages of Esfr {black squares) and model 
prediction (red triangles), for a set of spiral galaxies. For the model, H s = const, is adopted. 
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CO is not a linear tracer of the gas in GBCs, that £sf,gbc is not constant in these regions, 
or that age effects in the stellar population are impacting the estimate of S SFR . 



Blitz fc Rosolowskyl (120061 ) obtained an empirical fit relating the molecular-to-atomic 
gas mass fractions to a midplane pressure estimate, £ mo i/£ atom = (P^/P^o) 7 , f° r 



BR 



S(2Gp, 



(26) 



This pressure estimate assumes that the stellar disk dominates the vertical gravity, and 
combines the atomic and molecular gas into a single component. When combined with 
E SFR = £ mo i/£ SF , this yields 



J SFR 



g (Ph/Phfi 
t S F 1 + {Phi Pi 



h,0) 



(27) 



Blitz Rosolowskyl (120061 ) adopted a vertical velocity dispersion v g = 8 km s 1 , and ob- 



tained values 7 = 0.92 and Ph,o/k ~ 4 x 10 4 cm 3 K for the fitting constants. 



Leroy et al.l (120081 ) found that a similar relationship fit their sample of spirals, with 
7 = 0.8, Phfi/k = 1.7 x 10 4 cm -3 K, and a pressure estimate that includes gas self-gravity: 

ttGS 2 



Pt, = 



S(2Gp, 



'55 



the adopted value for the vertical velocity dispersion of the comp osite ISM is v 



(28) 
11 km s . 



Note t hat the grav i ty of d ark matter is not included in either the iBlitz fc Rosolowskyl (120061 ) 
or the iLeroy et al.l (120081 ) pressure estimate; dark matter may be increasingly important in 
far-outer disks or in low-surface-brightness disks. 

For Ph/Phfi 3> 1, equation (J27l) yields the same star formation rate as that in the GBC- 
dominated regime for the present model, Esfr — > E/t SF . For Ph/Ph,o 1, the limit in which 
atomic gas dominates molecular gas, equation (l2~7|) yields Ssfr S(P/j/p lj0 ) 7 /tsF- Since 7 
is close to unity for both of these empirical relations, and £ in the diffuse-dominated limit 
for both samples is typically ~ 5 — 10 M© pc~ 2 so that Ph varies approximately as pl^ 2 , 
these relationships are similar to the form of our result given in equation ( |22|) (see also eq. 
IA13[) . (Note that a value of 7 less than 1 partially compensates for varying S in P^.) As an 
example, Figure © presents the comparison between our model results and the empirical 
fits given above, for the galaxy NGC 5055. For this an d other gala x ies, w e find a close 
correspondence particularly with the empirical formula of ILeroy et al.l ( 120081 ) (as discussed 
above, slightly larger /„, /a than our fiduc i al cho ice shifts our predicted S SFR upward). The 
empirical formula of IBlitz fc Rosolowskyl ( 120061 ) produces somewhat more rapid decline in 
Ssfr at low S (for large radii) than the prediction of our model. 
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Fig. 7. — Same as in Fig. El for the galaxy NGC 5055 . Also included is t he comparison with 
the empirical formulae of iBlitz fc Rosolowsky ( 2006 ) (blue plusses) and Leroy et al. ( 2008 ) 
(green boxes), as described in the text. 
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4. Summary and Discussion 

1. Summary of the physical model - 

In this paper, we have developed a theory for self-regulated star formation in multiphase 
galactic ISM disks in which stellar heating mediates the feedback. The fundamental princi- 
ples we adopt are that for a time-and-space averaged steady state on ~ kpc scales, (1) force 
balance must be satisfied in the vertical direction, dynamically setting the midplane thermal 
pressure P t h of the diffuse gas based on the weight of overlying material (see eq. ITT]) ; (2) 
thermal equilibrium must be satisfied, with the heating rate set by the local star formation 
rate (see eqs. [13 [16]), and with the two-phase thermal equilibrium pressure -Ptwo-phase in the 
diffuse gas equal to the dynamically- imposed equilibrium pressure P t h; (3) the star formation 
rate is controlled by the amount of gas in GBCs (see eq. [[9]), with the (complementary) 
non-self-gravitating amount regulated by the thermal pressure (see eq. [H]). 

The set of algebraic equations embodying the principles above can easily be solved 
numerically to obtain Esfr as a function of the total gas surface density S and the midplane 
stellar-plus-dark- matter density p S( j, as described in the Appendix. An approximate closed- 
form solution, representing the key result of this work, is given by equations (I22p and (I23p . 
In the diffuse-gas dominated regime, S SFR oc S^/p^d (see eq. [22J) - As a consequence, no 
single Schmidt-type relation S SFR oc S 1+p is expected to apply in outer disks, since y/fhd 
need not vary as S p . 

Physically, outer and inner disks are distinguished by which gas component dominates 
the mass - diffuse or self-gravitating. Where diffuse gas dominates, the mean pressure and 
density and hence the cooling rate are fixed by the weight of the ISM. The amount of self- 
gravitating, star-forming gas created is then tuned to provide the needed FUV heating to 
balance cooling at this density: S SF r oc P t h oc £ v /p^j. Where self-gravitating gas dominates 
so that Ssfr oc S, the specific heating rate is fixed. The cooling rate of the diffuse gas 
depends on its density, which is proportional to the surface density S diff (and to the vertical 
gravity); S^ifr must therefore adjust until the cooling rate matches the heating rate. The 
limited surface density observed for H I gas in the central regions of galaxies likely owes at 
least in part to the constraint imposed by matching heating with cooling in the diffuse ISM. 

2. Connection to previous work - 

Our theory makes use of some of the same concepts - such as ther mal and dynamic 



equilibrium - discussed in previous work, but with a different emphasis. iParravand (119881 ) 



suggested that star formation is self-regulated in such a way tha t the UV radiation it produces 



maintains P ma x,warm near the thermal pressure of the gas, and iParravano fc Mantilla! (Il99ll ) 



applied this model to the Milky Way by adopting a radial profile for the warm gas density. 
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Considering outer galaxies, lElme green fc Parravanol ( 11994 ) pointed out that star formation 
could be strongly suppressed in outer disks if the midplane p ressure f alls su fficiently far below 
-Pmin.coid that even locally- compressed regions cannot cool. ISchayd (120041 ). treating the UV 
intensity as a fixed parameter, suggested that star formation would have a threshold imposed 
by thermodynamics. Here, we propose that (low-level) star formation is able to extend out 
to large radii in galactic disks because P m in,coid adjusts to follow the outward decline in the 
midplane pressure, via a decrease in the UV flux (tracking the decline in the star formation 
rate). 

In our model, star formation is regulated such that UV radiation created by young stars 
heats the disk just as much as is needed for the thermal pressure in the diffuse gas to meet 
the requirements imposed by vertical force balance. The regulation process depends on mass 
exchange between self-gravitating and diffuse components of the ISM such that star formation 
at the required rate can take place in the bound clouds. Since the star formation rate in 
normal galaxies is proportional to the mass in gravitationally bound GMCs, we believe that 
this self-regulation mechanism is the physical basis for the rela tionship between mo l ecular 
surfa ce density and midplane pressure empirically identified by iBlitz &: Rosolowskyi (120041 . 
20061 ). Th e theoretical relationship w e obtain is in fact slightly different from the empirical 
formula of IBlitz &: Rosolowskyi ( 120061 ) (see also lLeroy et al.ll2008l ). They found that the ratio 
of molecular-to-atomic surface densities is approximately proportional to an estimate of the 
total midplane pressure for a composite ISM (see $3] for details), whereas here we argue 
that the surface density of gas contained in GBCs should be proportional to the midplane 
thermal pressure in the diffuse gas. While physically and mathematically different, the two 
relationships yield quantitatively similar values of Ssfr provided the atomic surface den sity 
is relatively unifo rm (as is true in the observations of IBlitz fc Rosolowskyi l2004i 120061 and 
Leroy et al.ll2008l ). the ratio of thermal to total pressure is relatively constant, and the gas in 
GBCs is mostly molecular. In particular, the empirical BR relationship is most sensitive in 
the atomic-dominated limit, where the predicted star formation rate (assuming constant-tsF 
in molecular gas) is similar in form to our theoretical outer-disk law, Ssfr oc Sy/p^j. 



3. Comparsion to observation, present and future - 

Initial comparisons of predicted star formation rates with observed values, based on 
azimuthally-averaged data for disk galaxies, show quite promising results. In particular, the 
predictions follow the observations very closely throughout the two large flocculent galaxies 
NGC 7331 and NGC 5055, out to 1.2r25- In these initial comparisons, we have not "tuned" 
the model parameters a, f w , or tgF a ^ all, but simply adopted the same values for all 
the galaxies in the sample. By adjusting the parameters, closer agreement between the 
model prediction and observations can be obtained in several cases, by shifting the overall 
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normalization of Ssfr- Adjusting the prescription for converting stellar surface density to 
volume density can also yield a closer match to the data. 

Current data sets exist that make it possible to extend the present comparisons in 
several ways, including using ~kpc resolution maps (also including local variations of the 
gas-to-dust ratio) rather than azimuthally-averaged data, and considering dwarf galaxies. 
For far-outer disks, the necessary averaging scale is likely to increase, due to the flaring of 
the disk. It will be interesting to test whether galaxies with strong spiral structure, when 
examined locally, are consistent with the steady-state theory developed here, or whether 
transient effects within spiral arms are too rapid for (quasi-)equilibrium to be attained. 

In the analysis of §2, we employ several parameters (e.g. f w , a, t SF , and the ratio of 
Solar neighborhood pressure to star formation rate), adopting fiducial values that are based 
on current observations and/or theoretical work. As more detailed ISM information becomes 
available from extragalactic observations (such as local values of the gaseous vertical velocity 
dispersion, and the proportions of atomic gas in warm and cold phases), it will be possible to 
assign observed values rather than adopted parameters as inputs for predicting S SFR within 
individual galaxies. With local measurements of the stellar vertical velocity dispersion in the 
outer parts of individual face-on galaxies, it will be possible to obtain more direct estimates 
of the stellar midplane density, rather than simply adopting a prescription for the stellar scale 
height to obtain p s from E s . From surveys of edge-on galaxies, it will also be possible to 
obtain accurate measurements of the correlations of stellar disk flaring with other properties, 
that could then be applied to more face-on systems statistically. With more detailed data 
sets, it will be possible to test consituent elements of the theory - including equations ([7]) 
and (1TB]) - as well as the overall prediction for Esfr- 

4- Opportunities for numerical modeling - 

It is of considerable interest to test the idealizations of this theory, as well as its pre- 
dictions, via detailed numerical simulations. Numerical models must have sufficiently fine 
spatial grids ( £j pc) that the vertical direction is well resolved, must include heating and 
cooling such that warm and cold phases are present, and must model energetic feedback 
(leading to both heating and turbulent driving) from star formation in self-gravitating clouds. 
Because turbulent mixing at extremely small scales (<C pc) c an affect /,„, it is importa nt 



to assess this effect using very high resolution simulations (e.g. lAudit &: Hennebelld 120051 ) . 



With simulations, the turbulent velocity disper sion and warm/cold fract ions in the 



diffuse medium can be self-consistently calculated (cf. iKoyama fc Ostrikerll2009al ). such that 
the dependence of the ISM "state" parameters a and f w on the (input) galactic environment 
"variables" (£, p s , pdm, Z) and (derived) star formation rate Esfr can be assessed. The 
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star formation timescale tsF in GBCs and the ratio of mean thermal pressure (or mean FUV 
intensity) to Ssfr are also in principle calculable theoretically, although the dependence on 
the fundamental environment variables may be fairly complex. 

While we have not discussed here exactly how the formation of GBCs takes place, this 
can in principle be strongly affected by the angular momentum of the disk, with disks having 
very low or high Toomre Q departing from ob served star formation and m olecular /atomic (or 
GBC/diffuse) relations at a given £ and p s ( Koyama fc Ostriker 2009aj |bh . Normal galaxies 
have had sufficient time to evolve (lowering £ and raising p s by converting gas into stars) that 
the fundamental dependence of Esfr on angular momentum and shear may, however, not be 
evident in practice. Thus, exploring a wide range of types of observed systems, and testing 
both realistic and unrealistic galaxy models with numerical simulations, will be important 
for revealing the processes that control star formation at the most fundamental level. 



5. Limitations and prospects - 

By adopting a fixed value of the star formation timescale tsF in self-gravitating clouds, 
the present model is limited to the regime in which star-forming clouds ha ye "normal" 
properties, similar to tho se observed in the disks of Local Group galaxies (e.g. ISheth et al. 
20081 ; iBolatto et al.l 120081 ). In particular, it is not applicable to galactic center regions or 
starbursts where S exceeds the typical surface density ~ 100 M Q pc 2 of individual mid-disk 
GMCs. In such high-S regions, molecular gas completely dominates the ISM, but because the 
average density is higher than in mid-disk GMCs - so that the gravitational time oc p -1 / 2 is 
shorter, the star formation rate per unit mass is expected to be higher than it is for "normal" 
GMCs (consistent with observations). In detail, the star formation rate per unit mass in 
bound clouds in high-E regions is also expected to depend on the turbulence level (which 
is h igher in starbursts) and whether star-forming clouds are collapsing or marginally bound 
(cf. iKrumholz fc McKeeHiooS iRosas-Guevara et alJl20oit IPadoan & Nordlundll2009f l. Thus, 
from a combination of these effects, S SFR /E GB c is not expected to be constant in high-S 
regions, and the present theory should not be applied there. 

We note that low-metallicity systems, where GBCs are less self-shielded and have larger 
atomic-to-molecular ratios than "normal" GMCs of the same size and mass, may also have 
values of £sfr/£gbc different fr om the constant va lue t<7p = (2 x 10 9 ) -1 yr we have adopted 
(based on the CO observations of lBigiel et al.ll2008l ). In bound clouds with a large atomic-gas 
proportion (and with the carbon mostly atomic rather than in CO), the mean temperature 
will be larger than that in primarily-molecular clouds. Because of the lower internal Mach 
number, turbulent compression would be less extreme than in colder, primarily-molecular, 
clouds, which could affect the fraction of gas that is able to collapse and make stars. Addi- 
tional observational and theoretical work is needed to evaluate how the star-forming efficien- 
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cies of bound clouds depend on relative amounts of cold-atomic vs. molecular gas, and also 
to explore whether the masses and/or total column densities of star- forming bound clouds 
in low-metallicity regions (far-outer disks of spirals, and dwarf galaxies) differ systematically 
from the properties of "normal" GMCs. Observationally, a difficulty (particu larly in low 
metallicity regions) is that significant gas can be "dark" (IGrenier et al.l 120051 ) - i.e. not 
observable either in HI (because the hydrogen is molecular) or in CO (because the carbon is 
atomic). For Sola r metallicity and clo uds with A v ~ 8, the fraction of "dark" gas is expected 
to be only ~ 0.3 jWolfire et aDboiol ). 



The present model does not address radiative transfer effects explicitly. In particular, 
we have assumed that the optical depth through the diffuse gas is modest, such that the 
mean UV intensity and therefore the thermal pressure is approximately proportional to 
the local Esfr (see eqs. [TS] - ITS]) . Although this approximation becomes invalid where 
Sdiff ~ 20/ Z' M pc~ 2 (for Z' the metallicity relative to Solar), typically E 3> Ediff by this 
point (see e.g. Fig. [2]), so that changes in the FUV intensity and hence -Ptwo-phase and P t h 
would not significantly affect the predicted Esfr- An accurate determination of Edig where 
HZ' is high would, however, require an explicit radiative transfer calculation to assess G' for 
a given value of Esfr- This would be affected both by the amount and vertical distribution 
of diffuse gas, and by radiative transfer within star-forming clouds themselves. Radiative 
transfer effects also become important in the far outer regions of galaxies, where heating 
from nonlocal UV can exceed the local contribution. Results of radiative transfer models 
could in principle be tested by comparison to multiwavelength IR obs ervations, sin ce the dust 
temperature is sensitive to G' (i.e., to the mean UV intensity Jfuv) flDraind 120071 ). Here, we 
have not attempted to address these issues, but instead we have simply adopted an empirical 
Solar-neighborhood value for the ratio of Jfuv to Esfr to calibrate our relationships. 

Determination of the relative proportions of atomic and molecular (and, for low E, 
io nized) gas also d e pends on r adiative transfer. For GBCs , the solution for spherical clouds 
of iKrumholz et al.l ( 12009a ) and lMcKee &: Krumholzl (2010|) predicts the molecular-to-atomic 
ratio is (M atom /M mol ) cloud « [^ /a8 (^H,cioud/l-8 x 10 21 cm' 2 ) - 0.7]- 1 ; if clouds either have 
fixed total hydrogen column density or there is a known relationship between the mean 
surface density of this gas component S GB c (averaged over ~ kpc scales) and the column 
density of individual clouds, then this could be used to compute the contributions to E mol 
a nd Eat.nm from Eqr c . Th e diffuse gas could also be partly molecular; using the results of 



McKee fc Krumholzl (120101 ) for a slab of cold gas illuminated on both sides, a layer begins 
to become molecular when E co i d 11 M pc~ 2 /[Z'°' 8 ]. 

Finally, we note that simple models of the kind we have developed here - if validated 
by detailed numerical simulations and confirmed by observations - potentially provide a 
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valuable tool for studies of galaxy evolution. We caution, however, that careful appraisal 
of metallicity effects will be required before applying this model (or a refined version) to 
star formation in galaxies at high redshift. Our results could also potentially be adapted 
to provide subgrid ISM/star formation prescriptions for use in cosmological simulations of 
galaxy formation, but this too should be approached with care since the ratio of diffuse 
to gravitationally-bound gas is computed not for a local three-dimensional zone but for 
a vertically-integrated disk. Even with a subgrid model, the scale height of the diffuse 
warm+cold ISM (H ~ cr 2 z /g z ~ CTz/V^Gptot) must be resolved by several zones in order 
to obtain an accurate estimate for the midplane pressure, which controls the amount of 
star-forming gas. Because formation of gravitationally-bound star-forming clouds depends 
on non-local dynamical processes, simulations that are either resolved sufficiently to capture 
gravitationally-induced vertical motions, or are completely vertically-unresolved with a suit- 
able prescription for balance of ISM components, can represent the relevant physics more 
faithfully than simulations which resolve the disk with just a few zones. 
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A. Appendix 

Here, we provide an explicit formula that is solved numerically to obtain the fractions of 
diffuse and gravitationally-bound gas. We show how this formula leads to the equations (122]) 
and ( |23|) given in the text for the star formation rate. We also discuss how the galactic envi- 
ronment variables (£, p s d, Z) and the ISM model parameters control the transition between 
the diffuse-dominated and GBC-dominated regimes, and how presence of both regimes in a 
given annulus (due to spiral structure) affects estimates of the star formation rate. Finally, 
we provide a more stringent upper limit on Ejiff than the limit given by equation ( I24p in the 
text. 
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We start with the diffuse-gas thermal equilibrium equation (1X8]) relating the thermal 
pressure to the star formation rate E SF r = EgbcAsf> which may be expressed as 



1 Pi 



th 



th,0 



J GBC 



(Al) 



0d ^>SFR,0 tsF 

Here, P t h,o an d ^sfr,o are the Solar-neighborhood thermal pressure of diffuse gas and star 
formation rate that we adopt for our normalizations, and 

0.4" 

,' /. , • , \ 

1 + 3 



>d = T 



So 



(A2) 



is defined such that it is equal to unity for E = E , the Solar neighborhood diffuse-gas 
surface densi ty. The ratio of P t h,n to the mean local FUV intensity Jfuv.o is computed 
theoretically (jWolfire et al.l 120031 ). so that Solar neighborhood observations effectively enter 
the model through the ratio Jfuv,o/Ssfr,o- Using Solar neighborhood values Pth,o/k = 
3000 K cm and S S fr,o = 2.5 x 10" 9 M Q pc~ 2 yr" 1 with t SF = 2 x 10 9 yr as determined 
from extragalactic studies, this yields Pth/k = 600 K cm~ 3 0^" 1 (EGBc/ M pc~ 2 ) for the 
relation between thermal pressure and the surface density of gas in gravitationally-bound 
(star-forming) clouds]^ 



Next, we define 



E diff = xE, (A3) 
E GBC = (l-x)E, (A4) 

for the diffuse-gas and GBC surface densities, substituting these expressions together with 
equation ( lAlj) for P t h in equation fl9]) to obtain: 
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8nG( d clf w ap sd 



1/2- 



(A5) 



Here, we have introduced 



2aR 



th,0 



(A6) 



as indicated by Solar neighborhood observations of molecula r 



would yield a pressure a factor of ~ 2 below the observed local value. This simply reflects the fact that 



6 Note that inserting 
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Bronfman et al. 
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2 M P 



a_£C_ 



Zsfr.,o Z 10- y M per 2 yr- 1 if £ GBC 

of Ssf_r,o are higher by a factor ~ 2 dBertelli fc Nasil 12001 



a nd tsf = 2 x 10 9 yr, whereas observational est imates 



Vergelv et al.1 2002; F uchs et al 



2009). 



It is 

uncertain whether this discrepancy is the result of an underestimate of the local SgbCi an overestimate of 
the local EsfRj or a lower local value for tgp than the mean extragalactic value. 
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91 M pc" 



a 



Pth,o/k 



3000 K cm 



-3 



2.5 x 10- 9 M Q pc" 2 yr" 1 



^SF 



2 x 10 9 yr / 



Since -Pth,o/(0d£sFR,o£sF) = -Pth/ScBC and aP t h is the total (effective) midplane pressure, 
is the value that S would have to attain in order for the total pressure to equal 7tG££gbc/2; 
that is, for the pressure to be dominated by the gravity of the gas. 



Next, we define 
S = 
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(A7) 
(A8) 



which measures the relative importance of external gravity (stellar-disk + dark matter) to 
gas self-gravity in setting the diffuse-gas pressure. We note that if p st } is dominated by stars 
with p s = K/(2H S ) = txGY?J(2v\ 



J 2 ) 

J Z,SJ 5 



s 



(A9) 



which is oc (Q 9 „, s / 'Q s ) 2 in t erms of the stellar and gas Toomre parameters (cf. equation 4 in 
Koyama fc Ostrikerll2009bl ). Using fiducial values a = 5 and f w = 0.5, S ~ 16 in the Solar 



neighborhood. 

We can now re-express equation 



5J) in terms of the dimensionless variables S and 
ir = — (A10) 



yielding 
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Given values of p S( j and E in a galaxy, the variables S and w are set, and we can solve the 
(cubic) equation (lAlip for x numerically. The root x is bounded by and 1, so we use the 
bisection method. Given a solution for the diffuse-gas fraction x, the star formation rate is 
then 



J SFR, 



:i-x) 



(A12) 



We use this numerical solution to obtain the star formation rates for both our idealized galaxy 
(Fig. |2]), and for the comparison of the model prediction to the observed star formation rates 

(Figs, i- ED. 
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One important limit is that in which the ISM is dominated by diffuse gas, in which case 
x — > 1, and equation flAlip yields 

1 _ x _> w l + [l + ^] 12 ^ + 5l/2) _ (A13) 

When multiplied by E and divided by t$p, this yields equation (122]) of the text. 

Equation (I A 1 3 P requires w(l + S 1 ^ 2 ) <C 1 for self-consistency, whereas x — > in equation 
( lAllj) requires w(l + S* 1 / 2 ) ^> 1. An approximate solution for 1 — x allowing for both limits 

iS 1 1 

1 — x w(l + & ' ) 

The inverse of this, when multiplied by E and divided by t$F, yields the approximation for 
the star formation rate given by equation fl23l) of the text. Since w ^ 1 for the regions we 
are considering in this paper, it is generally the value of wS 1 ^ 2 (which depends on p s d but 
not on E) that determines which star formation regime holds. The approximation for 1 — x 
given in equation (IA14I) is good to within 11% for S > 10 and w > 0.01; for w < 0.2 
and 0.1 < S < 10, a better approximation (good to within 12%) is obtained by using 
0.5[1 + (1 + AS) 1 ' 2 } instead of 1 + S 1 ' 2 (cf. eq. lAl3il . 

Using equation (1A14D and the definitions of w and S, the overall timescale for gas to be 
converted to stars is given by 

^con = t; — (A15) 

^SFR 1 — X 



D / \l/2 
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dS SF R,o \2TrG( d clf w p sd 



+ t SF , (A16) 



where S 3> 1 is assumed for the latter expression. This timescale will be set by whichever 
of the two terms is larger. The first term is proportional to the vertical oscillation time, 
\/tt/ (Gp s d), which controls how fast cold cloudlets can sink to the midplane and form GBCs 
(for our fiducial parameter choices, this term in t con is 39 times the vertical oscillation time). 
The second term is the characteristic time for gas within GBCs to form stars. 

Quantitatively, the transition between the diffuse-dominated and GBC-dominated cases 
occurs where x — 1/2. From equation (IA110 . this yields the condition 

wy, = ^(1-Sw) = -\1- ^ j, (A17) 

where J2 h is given in equation (1A6[) . Since the right-hand side depends only p sd , this gives 
a value for the transition surface density E^ = Wi/2^h a s a function of p S( j. Taking the 
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fiducial parameter choices, Sw 2 = 0.37(p s d/0.1 M pc -3 ). For a given value of p S( j, the ISM 
will be diffuse-dominated if E < T,y 2 , and GBC-dominated if £ > Ex/2- For example, at 
the Solar circle where p S( j = 0.05 M pc -3 , the transition from diffuse-dominated to GBC- 
dominated would occur at £i/ 2 = 0.5411^ ~ 50 M Q pc -2 . We note that if p s d is large enough 
that Sw 2 > 1, then x < 1/2 (i.e. Egbc > Edm) regardless of the value of E. 

At a given galactocentric radius, if the surface density satisfies either £ <C £iy 2 or 
E ^> Ex/2 at all azimuthal angles, then a single star formation regime applies, and the 
azimuthally-averaged star formation rate can be obtained from the azimuthally-averaged 
gas surface density. If, however, there is a transition from E < £1/2 in interarm regions 
to E > Ex/2 in spiral arm regions, then the star formation regime changes from diffuse- 
dominated to GBC-dominated, and the prediction of E SF r based on (E) would depart from 
the true value due to nonlinearities. Equation (1A14[) . when multiplied by tsF and evaluated 
using the value of w in the arm gives t C on,arm = E arm /EsFR,arm in the arm gas: 

fs! ' +t SF (A18) 



' o nj .ann ' ^ -,1/2 



Warm "I" W arm >S'arm 

(note that 5 ^> 1 does not hold in arms, so that eq. IA16I should not be used; w aTm may 
however be £ 1). An analogous expression holds for t con ,ia = £ia/£sFR,ia in the interarm 
region using w aTm — > and S^m — > Si a . Letting / arm be the mass fraction in the arm in a 
given annulus, the star formation rate in the annulus using the arm and interarm conditions 
separately is 



J SFR,arm+ia 



^con,ia 



I 1 / ^con,ia -. 
~r /arm ( ~ J- 



'con, arm 



(A19) 



At a given radius, wS 1 ^ 2 oc p^ 2 varies by £ 10% from spiral perturbations, so that 



UnaS 1 ; 1 / 2 ~ w arm sUm — > wS 1 ^ 2 . For most regions of interest, w ia <C wS 1 ^ 2 C 1 so that 
^con,ia ~ tsp^wS 1 ^ 2 )^ 1 . Equation (IA19I) can be compared to the star formation rate that 
would be estimated using (E) in equation (IA14I) . which yields EsFR,az = (E)(l — x) az /tsF ~ 
(E)wS l/2 /t SF ~ (E)/t CO n,ia (assuming (w) < wS 1/2 ). Taking the ratio of E SFRiarm+ia to 
EsFR.az, we obtain 



' S "' ,{ - i,n " +la ~ 1 • /,„,; ( —-, f/ ' ar;: ' ,, ., , ) . (A20) 



E S FR,az ywS^lWarm + wS 1 / 2 + 1] 

Since the term in parentheses is typically order-unity, the true star formation rate (i.e. 
EsFR,arm+ia) can be considerably larger than the estimate E SFR az based on the annular az- 
imuthal average (E) if the gas is highly concentrated in the arms. 

Finally, we consider the upper limit on the diffuse-gas surface density. Taking the inverse 
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of equation (IA11I) and multiplying by E, we have for the diffuse-gas surface density 

(A21) 



^diff 



+ 



S(l-x) 1 (1-ai) 



y-2 



1/2 



< — T72- ( A22 ) 

l+(l + % + 5) 1/2 

The right-hand side of the inequality (1A22I) is the limiting value of Edifr for i C 1, i.e. 
the case in which the ISM is dominated by GBCs. An absolute upper limit Edis < S^/2 is 
obtained by taking E/j/E, S — » 0; the result is given by equation (1241) of the text. In practice, 
the terms E/j/E and 5* in the denominator of equation flA22j) are appreciable, so that Edis 
is below E/j/2 by a factor of a few. 
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